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Abstract 


A class of optimal linear time-invariant tracking systems, 
both in continuous time and discrete time, of which the mamber of 
inputs (which are restricted to be step functions) is equal to 
the naimber of system outputs, are dealt with in this report. Along 
with derivation of equations and design procedures, two discretiza- 
tion schemes are presented, constraining either the control or its 
time derivative, to be a constant over each sampling period. De- 
scriptions are given for the linearized model of the F-8C aircraft 
longitudinal dynamics, and the C* handling qualities criterion, 
which then serve as an illustration of the applications of these 
linear tracking designs. A siiboptimal reduced state design is also 
presented. Numerical results are given for both the continuous time 
and discrete time designs. 
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CHAPTER 1 
INTRODUCTION 

This report deals with a class pf optimal linear tracking systems, 
in continuous time and in discrete time. An application to 
the control of an aircraft is discussed, along with a suboptimal 
design. 

The type of designs that this report focuses attention on is 
•compensator design for a class of multi-input-multi-output linear time- 
invariant dynamical systems in which the nuxriber of controls is equal to 
the number of system outputs, and the system output vector is required 
to track a reference input vector of step functions with zero steady . 
state error. It is found that the resulting designs are conceptually 
the proportional-ihtegral-derivative (P-I-D) controller designs. 

Sandell [1] ,. and Sandell and Athans [2] , present a continuous time com- 
pensator design for such systems. Their designs suffer from the practi- 
cal drawback that the initial jump in the control is very often too fast 
and may rate-saturate the control actuator. We pursue this design in 
continuous time, but eliminate the above stated effect by imposing 
continuity of the control everywhere. Holley and Bryson [3] give a regu- 
lator design with non-zerb set points, and it is found to be of the same , 
structure as the one that we are ensuirig. Owing to the need 
for in 5 >lementing control systeina ;in 4^ we concentrate much 

of our effort into getting a parallel design of a discrete time compensator. 

. ' ’ ' - 8 - 



The NASA F-8C digital- fly-by-wire (DPBW) aircraft is a high per- 
formance aircraft and is considered a very convenient testbed vehicle for 
many' flight control systems. One of the purposes of this research is to 
investigate and to characterize the idea of designing a conceptual P-I-D 
controller for the linearized longitudinal dynamics of this aircraft, 
using a C* handling qualities criterion [10] , [11] . 

The C* parameter of an aircraft is a linear combination of the 
pitch rate and the normal acceleration (and possibly the pitch angle) of 
the aircraft. We shall discuss its significance and motivation later, 
along with a description of the linearized models of the F-8C aircraft 
longitudinal dynamics. Briefly we can view the design problem as a 
tracking problem [15] as follows: the pilot, by shifting the command stick 
position, inputs a desired steady state value of the C* parameter to the 
compensator, which is then required to track the system C* response to 
this value, by minimizing an appropriate penalty function, with zero 
steady state error. 

Athens, et al. [8] discvisses briefly the desire of using the 
C* criterion in the control of the F-8C aircraft. Honeywell staff [12], 
[13] discuss a C^-design using model following techniques. Unfortunately, 
their design procedures are not transparent enough to us, at least from 
a theoretical point of view. 

It is our desire that eventually our design procedure will be 
incorporated into the multiple model adaptive control (MMAC) scheme 
assumed in [8]* Henceforth, we follow a design philosophy similar to [8] 
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and use the same set of data that it employs. We ej^jloit this further 
in Chapters III and IV. 

The remaining chapters of this report are organized as follows. 

In Chapter II, we present most of the theoretical results required 
for the deterministic designs in this report, together with some dis- 
cussions on these theoretical issues. First we derive the equations 
for the P-I-D controller in continuous time, as a modification 
of the results in [1] or [2]. Then we obtain similar results in the 
discrete time case. Owing to physical motivations, we discuss two dis- 
cretization routines where we treat (1) the control, (2) the rate of the 
control, as a constant over each sampling period. In other words, method (1) 
uses a zero-order hold and method (2) a first-order hold. As we shall 
discuss later, due to sudden changes of the control, i.e. high rate of 
change, method (1) may not be physically satisfactory for aircraft imple- 
mentations. However, the difficulty of getting a proper linear realization 
of method (2) may also impose additional problems. 

After the completion of the work in this research, it has come 
to our attention that T.A.S.C. also came up to some p-I controller 
designs in [4] for VTOL aircraft control. The last part of this 
chapter is devoted to brief con^arisons of these results. 

In Chapter III, the general motivations of applying the C* 
criterion in the controller design of the F-8C aircraft linearized 
longitudinal dynamics are discussed. We present a review of the 
linearized models of the F-8C aircraft longitudinal dynamics (at 16 
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flight conditions) , and the C* criterion. We then formulate the C* 
tracking problem with reference to the results in Chapter II. 

Reference [8] presents a good approximation of the control gains 
and transient responses of the linearized F-8C aircraft full state 
longitudinal dynamics by a reduced state model. This approximation is 
also applicable to our present problem. The task of Chapter IV is to 
give an account of such a reduced state model. Theoretical comparisons 
between this model and the full state model are made. 

In Chapter V, we present some numerical results, both in con- 
tinuous time and discrete time. Deterministic designs, 
simulations, mismatching of flight conditions, and full state 
system versus reduced state system in performance are the major issues 
that we touch upon. 

Finally, in Chapter VI, we discuss the overall results and 
draw some final conclusions. Also we outline some possibilities for 
future research. 


CHAPTER II 


P-I-D CONTROLLER DESIGNS 

In this chapter, procedures are given for designing P-I-D con- 
trollers for multi-input multi-output linear time-invariant systems. 
Although results are first given for continuous time systems, the main 
parts of this chapter are devoted to discussions of discrete time 
systems, discretized from the continuous time systems via two different 
discretization schemes with sampling period T: 

(1) For each integer k, we constrain the control to be a 
constant over the saitpling interval [kT, (k+l)T) . 

(2) For each integer k, we constrain the time derivative 
of the control to be a constant over the sampling 
interval [kT, (k+l)T). 

For the discrete time systems obtained by scheme (1) , we derive 
three different sets of procedures for designing the P-I-D control 
gains. The first procedure is the discrete time equivalent of the 
continuous time design procedure; the other two are siitply approxima- 
tions of the control gains of the continuous time systems. 

One drawback in using scheme (1) is that it introduces a series of 
steps in the control, and hence, a series of impulses in the control 
rate, which might then lead to rate saturation of the control actuator 
in many physical systems. 

Scheme (2) is discussed in some detail. There, the derivation of 


the design procedures bears many similarities to that for scheme (1) , 
although required conditions are som-ewhat stronger. 

We divide tliis chapter into five sections. In Section 2,1, a 
design procedure for the continuous time case, which is almost identical 
to that given in [1], is presented. Sections 2.2 and 2.3 are discussions 
on the discretization scheme (1) , the first of which deals with the 
scheme itself and a direct design procedure of the discrete time 
systems, while the latter is concerned with approximations of the con- 
tinuous time design. Section 2.4 deals with discretization scheme (2) 
and Section 2.5 contains some final remarks and a coitpariscn with the 
results in [4] . 

2 . 1 Continuous Time Designs 

■ ■ ' '■ * I 

As a first result, we present a design procedux'e in the continuoiis 
time case, which can be proved by adding in a constraint that the 
control has to be continuous for all time, and then using the method 
employed by [1] or [2] . For this reason, the proof is omitted here. 

A derivation of this procedure, which is of a somewhat different 
flavor, can also be obtained if we attack this problem by using the 
results in [3]. 
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Proposition 2.1 Suppose that we are given a linear time-invariant 
dynamical system in continuous time: 


x(ti) = A x(t) 
v{t) = G x(t) 


+ B u(t) I 


System I 


( 2 . 1 ) 

( 2 . 2 ) 


, . n . . ■ m , , m 

where 3c(t) e R , u(t) G R , ;^(t) e R , and A, B , C are n><n, nxm, mXn 
real constant matrices respectively. 

We assume that 

(i) rank A = n 

(ii) [A, B, C^] is a minimal realization 

(iii) rank £ ^ = m. 

We wish to minimize a cost functional; 

I oo ... ... . 

J(u) = y [(Z^ - X.^T) ) ' ) + u' (T)R u(T)ldT (2.3) 

0*' I ^ 

where z is an arbitrary step input, and Q > 0, R > 0 are rP<n 

~o ■ ' , . — 

positive semi-definite and mxm positive definite matrices respectively. 
Then there exists a control u(t) such that the , performance index J 
is minimized. The inputs are tracked with zero steady state error 
for all step inputs, and the closed loop system is stable. The con^- 
trpl is given by 


u(t) = j ~ 


- X ) + u(0) 
-o — 


(2.4) 
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where 


L = (K - K a"^B) (C 

M X 


N = (K^ + it £.) A 


-1 


(2.5) 

( 2 . 6 ) 


The matrices are found from the positive definite solution of 

the algebraic matrix Riccati equation 


P P 

^1 i -12 


P' P 

-12 -22 


P P 

“11 ^12 


pi p 

-12 -22 


A B 


0 0 


A' 0 


B' 0 


P P ' 

-11 -12 


P 't P 

-12 -22 


r"’''[0 I] 


£ll 


p I 

LH2 


(2.7) 


^12 


p 

-22 


£'R £ £l 


0- 


with 


[K^ [0. II 


p 

hi 1 

1 ^12 

1 

p ' 

,-12 

-22 


( 2 . 8 ) 


where P , P_ are matrices of dimensions n><n, nxm/ mxm respectively 

.XX Xa 

The structure of such a control system is shown, in Figure 2,1, 




Fig. 2.1 • Continuous Time Linear Tracking (P-l-D) Compensator Design 



Note that in equation (2.4) ^(0) can be arbitrary. In particular, 
there are the following choices. 

The choice of u(0'*’) in [1] optimizes the initial value to the integrator, 
in the sense that the cosr'fundtional in equation (2.3) is minimized with 
this particular choice. However, this procedure does not pay any regard to 
u(0 ), or to be specific; ’ the rate of change of the control at time t = 0, 
which is not penalized in the performance index at all. Such an initial 
jump in control may rate saturate the control system and lead to 
subsequent • failure of the future control * On the other hand, [2] 
picks a choice such that the control, u(t) , at all time t is independent of 
3c(0) . The problems of rate saturation are also present here, again due to 
the difference between u(O^) and u(0 ), In order to eliminate rate satura- 
tion, we need to impose Continuity Of the control, in particular/ if we 
constrain 

uCo"^) = h(0“) : (2,9) 

Then it is easy to see that the control is continuous for all time t. This 
is indeed our choice for this report . 

Now we present a similar design procedure for a discrete time 
system with a similar performance index . 

2,2 Discrete Time Designs, scheme (1) 

In Section 2.1, a design procedure for continuous time systems 
is discussed. A similar procedure can be derived for a class of dis- 
crete time systems, namely, systems with a constant control over each 
sampling period, and a Similar performance index. We discuss this in 


I 
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the following proposition. Then we surnmarize the relationship between this 
result and the previous proposition. 


Proposition 2.2 Suppose that we are given a linear time-invariant 
dynamical system in discrete time (with sampling period T) 


x^Oc+DT^ = 


;^(kT) = ^x(kT) 


( 2 . 10 ) 


( 2 . 11 ) 


where x(kT) e u(kT) e r”', 


e R , and A, , B, , C . are n^n, nxm, 
'—a —a 


mXn real constant matrices respectively. 


We assume that 


(iv) rank (A, - I) = n 

(v) [A , C ] is a minimal realization 

“Ti ~Q ~~a 

(vi) rank ^ (^ - I) ^B, = m. 

-ni — d — • — d 


We want to minimize a cost functional, 


V- ' X.tiT))'2^(z^ - ii(jT)) 


+ (U( ( j+1) T) - U( jT) ) R^ (u( (j+l) T) 


( 2 . 12 ) 


- u(jT))j ‘ 


where z^ is an arbitrary step input, and 2^ ^ > 0 are nxn 

positive semi-definite cund mXm positive definite matrices respectively. 


Then there exists a control sequence such that the performance 
index is minimized. The inputs are tracked with zero steady state 
error for all step inputs and the closed loop system is stable. The 
control is given by 
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with 



where » £22 matrices of dimensions nxn, nXm, mxm respec- 

d d d 

tively. _ , 

Proof 

^ Let ■■ . ■■■ , 

x(kT) = x(kT) - X (2.18a) 

— ■ ■ ■ . ■ ' , ''-—O 

u(kT) = u(kT) - u (2.18b) 

— ■ — ■ ' ■ 

v(kT) * u((k+l)T) - u(kT) (2.18c) 


where 


X - - 


u 

“O 


(Sa - P '^SaHo 
<Ca<^ - £)■%> '"Sa 


(2.19) 

( 2 . 20 ) 


which is easily seen to be the unique solutibn of 
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X = A, X + 1^ U 
— o ~d — o ~d — o 


z = c,x 
— o — d— o 


(2.21a) 

(2.21b) 


Then the dynamics of the system can be written as 


x( (k+l)T) 




! 

■ i 

x(kT) 

+ 

'o_' 

u ( (k+1) T) 


g_ 



u(kT) 

m “ 




V (kT) 


( 2 . 22 ) 


and the performance index given in equation (2,12) is equivalent to 


J = ^ |[x* (jT) u' (j 
j=0 




io 


+ 2L* UT) ^ 


X(jT) 

u(jT) 


(2.23) 


Hence the problem is reduced to an optimal linear quadratic regulator 
problem, since > 0,^ > 0, by [ 5]» there exists a solution 

to this regulator problem given by 


v(kT) = (K K, 1 
d '^d 


x(kT) 

: I 

u(kT) 


(2.24) 


where are obtained by equations (2.16) and (2.17) . 

a ^d 
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Note that 


u(kT) = x((k+l)T) - Pi, x(kT) 
—a — — . — 


(2.25) 


and assumption (vi) implies that B, is of .rank m. Therefore there exists 

-d . . 

a mxn matrix F, such that 
—a 




(2.26) 


which gives us 


u (kT) = X ( (k+1) T) - X (kT) 


(2,27) 


Substituting this into equation (2.24), we get 


v(kT) = (K, - K + K* F, x((k+l)T) (2.28) 

- 1 ^ - 2 ^ -^-d - -2^-d - 


or equivalently, 


V (kT) = (Kj^ - K 2 (^ - 1) ) X (kT) + ^2 (k+1) T) - x(kT) ) 

(2.29) 


Now choose Lj such that 
.."+4 ■ 


£1 - & i:a<5a - i> * -ia£d 

d:, a 


(2.30) 


which gives 
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(2.31) 


But from equation (2.26) 




(2.32) 


Therefore , 


■ id = % - ill <Sd - i>'\> <£d‘2d - 

a d 

Using equation (2.30) and noting that 


^ “ 5 , 


-o 


(2.34a) 


and 


x((k+l)T) - x(kT) = x((k+l)T) - x(kT) 


(2.34b) 


we get 


x.e. 


v(kT) = Ljj(z^ - X.^kT) ) + K2 (x( (k+l)T) - x (kT) ) 

d ■ 


u((k+l)T) - u(kT) = u (z - ;^(kT) ) + k (x( (k+1) T) 

^ 'Cl 


(2.35) 


- x(kT) ) 


(2.36) 


-24- 


which can also be v^ritten as 
k-1 


u(kT) (z^ - ^(jT)) + ^(x(kT) -x(0)) + u(0) (2.37) 

j=0 


where 


d 


(2.38) 


and this completes the proof of this proposition. 

Notice how similar this solution is to the one in proposition 2.1. 
The structure of the compensator thus constructed is shown in Figure 2.2. 

A natural question would be what equivalence can one get between 
this discrete time system and system I . We answer this question in the 
following proposition. 

Proposition 2.3 Suppose that we discretize system I in proposition 2.1 
with sampling period T such that for each integer k, the control u(t) is 
constrained to be a constant over the sampling interval [kT, (k+l)T), i.e 
for each t in this interval/ 

u(t) = u (kT) ; kT£t£(k+l)T (2.39) 

Then we have a sample-data realization as follows, 

I 

ir( (k+1) T) = ^x (kT) + a(kT) A i (2.40) 

/ System II 

• ;^(kT) = x(kT) - (2.41) 




25 - 
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with ^ matrices given by the relations 

-nQ — O 

(2.42) 

(2.43) 

(2.44) 

where A, C are the matrices defined by system I, eqs. (2.1),(2.2) . 

If in addition, we have for all integer k 

ImlX^(A) - X^(A)]?^k^ (2.45) 

whenever 

Re[X^(A) “ Xj(A)] = 0 (2.46) 

for i y j , where (A) , X^ (A) are (distinct) eigenvalues of A, then 
the assumptions (i) , (ii) , (iii) , in proposition 2.1 imply the assump- 
tions (iv) , (v) , (vi) in proposition 2. 2. 

Proof 

Here we only prove the validity of assumptions (iv) and (vi) due 
to (i) and (iii). The other parts of this proposition are all contained 
in [6 J, hence their proofs are omitted. 

Since assumption (i) , i.e. 

Rank A = n 


, AT 


2d 




,-l at 

= A (e— - I) B 
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implies that for all i, 0 < i 5. 


X^{AT) 0 


therefore 


. 1^0 


But for every 0 < j £ n, there exists an 0 < i £ n such that 

X.(AT) 

■V / AT ■ . i'— ' . 

X . (e— - I) = e - 1 


x.e. 


rank (A, - I) = n 
— d — 


which proves assumption (iv) . 
Now note that 


-*1 AT —1 —1 AT 

‘ A ^ (e-^ -1)B. 


= C a“^b 


Assumption (vi) is then obvious from (iiLi) , 

From now on we shall refer to this discrete time system as system II 
Note that this system does not give us an equivalentt performance index 
(equation (2.12)) to that in system I (equation (2.3)) . 

t equivalent to the extent that iu(t) is kept to be a constant over each 
sampling interval while evaluating the cost functional in equation (2.3) 
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2.3 Two Approximate Schemes for Discretizing the Continuous Time 
Compensator Design 

The objective of this section is to investigate two ways of how 
the continuous time linear .tracking compensator design presented in 
Section 2.1 can be discretized into discrete time designs by making some 
approximations. 


The discretized system that we want to characterize here is one 
which has a constant control over each sampling interval. To be specific, 
given a continuous time system and a control law which tracks the system 
output to a step input, z^. 


x(t) = A x(t) + B u(t) 

(2.47) 

^(t) = C X(t) 

(2.48) 

/•t 


u(t) = L 7 (Zq - ^(T))dT + N(x(t) - x(0)) + u(0) 

(2.49) 


where L, N are matrices specified by equations (2. 5) -(2. 8), we propose 


to find a control law with control gains , given by 

k-1 

u(kT) = L 2 - u(jT)) + N (x(kT) - x(0)) + u(0) 

-D ^ " 

u(t) = u(kT) for kT < t < (k+l)T 


(2.50) 

(2.51) 


which I represents an approximation of the continuous time control law 
(equation (2.49) ), and is applicable to the sampled dynamics of the 
system, with sampling period T, 
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AT 

3c ( (k+1) T) = e— x^(kT) + 
bo 


e— ^dx B u(kT) 



(2.52) 


^(kT) = C x(kT) 


(2.53) 


We proceed to find ^ by two different sets of approximations: 

(i) First-order approximation 

We assvime that approximately. 


J rkT 

0 


k-1 

% - J 1 (t) )dx « T L 2 (Zq - Z( 3 T) ) 


(2.54) 


Here as a first-order approximation, we take 


L = T L 
-D — 

N - N 
-D — 


( 2 . 55) 
(2.56) 


(ii) Approximation with constant integral over sampling interval and 
pole-allocation 

Noting that R, and are feedback gains on the states in the 
continuous time and discrete time systems, respectively, we try to make 
use of the pole allocation method. In closed-loop form, the system 
dynamics can be written as 


X(t) = (A + B^^^N)x(t) + 


BL r 

Jn 


~ ^ ^ + B(u(0) -N X(0) ) 


(2.57) 


for the continuous time, and 
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k-1 


X(k+1)T) = + Bj^)x(kT) + E (5o ■ L(:iT)) 

j=o 


+ B_(u(0) - Nx(0)) 

"■jj '** 


(2.58) 


for the ^discrete time. We assume that 

Jt ‘ , 

j (z^ - X.(T))dT w T E 


(2.59) 


for kT < t < (k+l)T, and 


B(u(0) - N x(0)) = B (u(0) - Nx(0)) 
-D — -D— 


(2.60) 


which is exactly true when u(0) and x(0) are zero. Therefore, we may 
now write 


* a. o XT (A + B N)T 

Ad + S£&=®'” ““ 


B L 
-D-d 


•'0 


(2.61) 

(2.62) 


which give us, since B, hence is of maximal rank, a (minimum norm) solution 






"D =-D 


' •'o; 


(2.63) 

(2.64) 
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2.4 A Discretization Scheme Which Reduces Rate Saturation 

In the last two sections, we discussed a discretization scheme 
which constrains the control to be a constant over each sampling interval, 
and some compensator designs associated with it. These designs, however, 
are susceptible to control rate saturation because the derived controls 
are step fvmctions. As an altemati^>e, in this section, we present a 
second discretization scheme which constrains the time derivative of the 
control to be a constant over each sampling interval, but enjoys an 
equivalent'^ performance index to that in system I, subject to this con- 
straint. This has the merit that the problem of control rate saturation 
can be reduced considerably. 

Proposition 2 . 4 If we discretize system I with scunpling period T by 
constraining the time rate of the control, u(t) , to be a constant 
over the sampling period [kT, (k+l)T) for each integer k, i.e. 

u(t) = u(kT) = constant for kT £ t < (k+l)T (2.65) 


then we get an equivalent realization in discrete time, 



'^'equivalent to the extent that u(t) is kept to be constant over each 
sampling interval, while eyaluating the cost functional in equation (2# 3) . 
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where ^ are matrices from system II, given in propositions 2.2 
and 2.3, the £ matrix is from system I , and is given by 


T 


=y y* e-’’^dT ds B 


( 2 . 68 ) 


o o 

where A, B are from system I. 


I - 


Moreover, we obtain an equivalent performance index,' 


J = ^ |[(x(jT) - ^)' (u(jT) - u^)'] £ 
j=0 


x(jT) - X 
— — o 


u(jT) - u 
l_— ~o J 


+z[(x(jT) - - u ) ]M u(jT) 


(2.69) 


where 


+ u( jT)R u(jT) 


with 




^ At,^ 
i— C'2^ C er- dt 


M = 


t s 



(2.70) 


(2.71) 


R = R T 


{2.72) 


" -(Sa - i>'^2dHo 


^ = -(C(A - I)“V)"^Z 
-o — -d — -d — o 


(2.73) 

(2.74) 
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which are the steady state values of x(t) and u(t) in both continuous 
time and discrete time. 

Proof 

Firstly note that ^ given by equations (2.73) and (2.74) 
are unique solutions to any one set of the three following sets of 
equations, thus estciblishing the equivalence of steady state values 
in continuous time and discrete time. 




0 A X + B u 

— — . — o — — o 


z = C X 

— o I — — o 



A,x + B,u 
-^-o -d-D 


C 


X 

-o 


(2.75) 


(2.76) 



\ 



X 

-o 

+ 

D 

— o 

0 

I 


u 


T I 


■ 


”0 




(2.77) 


0 


It is easy to see that equations ( 2 . 73 ) and (2.74) indeed give the 
unique solution of 3, and 3 and y have the same solutions for x^,u^. 
Therefore, all we need to do is to show that they are also solutions 
to a. Solving a, we get 


X = -A B u 
-o — — —o 


u^ = -(C a“^b)”V 

“O — — — — o 


(2.78) 

(-2.79) 


-34- 


But from equations (2.42) and (2.43)# we know that 


(^ - 


(2.80) 


which implies that equations (2.73) and (2.74) are also solutions 
to a. 

System I can be realized as 1 


A. 

dt 


x(t) 

u(t) 


X 

—o 


u 

-o 


B 


^(t) - z = Ic a] 


JCu) 


=/{[U 


x^(t) 

- X 


”Q 

u(t) 

- U 


— o 



x(t) - 

X 


~o 

u(t) - 

u 


-o 


r 


u(t) (2.81) 


(2.82) 


[{x(t) - X ^ (u(t) - u.) 1 
— — o — — o 


£'2. £ £* 


x(t) - X " 
— “O 

0 a 


u(t) - u 

- ■ 


•“ “O 


+ U' (t) R Ji(t) I dt 


(2.83) 


Constraining u(t) = u(kT) for kT £ t < (k+l)T, we get by [ 7 ] an 
equivalent discrete time realization# 
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x( (k+l)T) - X 


u( (k+l)T) - u 
“ ~~o 


= exp 


A B 


0 0 


x(kT) X 
— -o 


u(kT) - u 
— _o 


+ / exp 


A B 


0 0 


u(kT) 


(2.84) 


jr(kT) - i = [C 0] r x(kT) - X 

“■ “ I — — o 


u(kT) - u 
— 


(2.85) 


and an equivalent perfoannance index 


■SlH 


(x(jT) - x^) (u( jT) 


u^) 2r2L<3T) - ^ 


Ut^T) - U 
— -O 


+2[(x(jT) - X ) (ui 


\^) ]M U( jT) 


( 2 . 86 ) 


+ u(jT) R u(jT) 


T (\ ±' 


£'2.G 0 


0 I 0 


A B 


0 0 


(2.87) 


£'£ £ £ t / A B 


0 0 


[0 0 j ■ Li 

( 2 . 88 ) 
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R = RT + [0 I] 

o o 


X [ J exp 
o 


/ 


t 

exp ■ 


”a' 

o' 

sj 

1 

^ds 

C'£ C 

£ 

B' 

0 


\ 

0_ 

0 






- 


Fa 

b1 


2. 


s>ds]dt 


0 

0 

) 

I 

L~ -J 




After some algebra, we see tlat 


exp 


'a 

-1 

J . 
* ( 


2d‘ 

0 

£ 

j 

0 

£ 

- 

J 





I 

/ 


T 

exp 


'a 3 


'o' 


D 

— 

) 



— o 


T|dT 


= 


0_ 1 

) 

L 


T £ 

•m - 






(2.89) 


(2.90) 


(2.91) 


£ = 2 


M = M 


(2.92) 

(2.93) 


R = R (2.94) 

' A- 

where M, R are given b_T equations (2.68) , (2.70) -(2.72) , thus 

completing the proof of this proposition. 

Rrom now on, we shall caJ, this system, given in proposition 2,4, 
system III. We are now in a position to obtain a compensator design 
such that the performance indsx J in equation (2.68) is minimized .s. It, 
tvurns out that we get a very similar expression for ii(kT) as in propo- 

■ " ^ ' - ' , . ' 'i- - 
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sition 2.2 for system II. The results are given in proposition 2.5. 
Since the technique in proving them is similar to the one we used 
in proposition 2.2, we only give the major steps in the proof.’ 





-38- 


b) If in addition, we assume that 


(viii) rank [A, + D K, - I] = n 

— 0—1 — 


- 1 , 


(ix). rank [C(^ ^ 5 d ^ 2 ^ ^ * 


m 


then the control u(kT) is given by 
.k -1 ' * 

u(kT) = ^ ^ " X.^ jT)) + N(x(kT) - x(0) ) + ij(0) 

. j=o 


(2,98) 


where 


- 1 , 




xlC(A + B K - I)"^(B. + D K.)!*'’ 
d ^o-i — -d 


(2.99) 


N = (T K, + L C) (A, + D k, I) 

— —1 — — — O-i — 


-1 


( 2 . 100 ) 


Proof 

For (a), the result is immediate by using [ 7l* For (b) , from 
equation ( 2 . 66 ) and (2.95), we obtain 


(x((k+l)T) - X ) = [^ + D K, Kx(kT) - X ) 
— — o -^d —0—1 — — o 


t [B, + D K^] (u(kT) - u ) 
-d - 0—2 — -o 


( 2 . 101 ) 


Assxiraption (ix) implies that there exists a mxn matrix F such that 


f(b^ + 0 ^ 2 ) “ L 


( 2 . 102 ) 
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Hence using equations (2.101), (2.102) and (2.66) in equation (2.95) , 
we obtain 

^ (U((k+1)T) - u(kT)) = [K- - K-F(^ + D K-)](x(kT) - X ) 

T — — ■ —1 —2 ^ —0—1 — — O 

+ K_F(x((k+l)T) - X ) (2.103) 

—2— — — O 

This implies 

k-l‘ / 

u(kT) = t[K^ - Soil - I) ^ - ^) 

• j=0 

+ T K 2 F(x(kT) - x(0) ) + u(0) (2.104) 

After some manipulations, this reduces to equation (2.98) and thus 
completes the proof. 

Note that in synthesis, this compensator design has the scune 
structure as the one for system II . Also note that the validity of 
assumptions (viii) and (ix) , unlike assumptions (i) - (vii) depends 
hot only on the system matrices, but also on the weighting matrices in 
the performance index. Furthermore, equation (2.66) gives 

X({k+1)T) = ^x(kT) + (B^ - i- D^)u(kT) + |- D^C (2.105) 

The presence of the last term on the right side of equation (2.105) appears 
to be a serious drawback in using the control law which we have 
just formulated in proposition 2.5, because computation of ^( (k+l)T) in 
this scheme requires knowledge of jc ( (k+1) T) . It appears that the system has 
become noncausal , and some approximations must be used to make it causal. 
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TWO possible ones, which are both ad hoc, are 

a) approximate u(kT) by ^ (u (kT) ^ u(k-l) T) ) , giving 

X((k+1)T) = A x(kT) + (B + u(kT) - U((k-1>T) 

(2.106) 

If T is large, the deviation introduced by such an approximation might 
be considerable, especially for the first few time steps, where 
the system is by no means near steady state. 

b) sinply use 

x((k+l)T) = ^x(kT) + ^u(kT) (2.107) 

— ^ —d— -d— 

but again, this may lead' to rate-saturation. 

Depending on the actual system, any such approximation 
can induce instability, a highly undesirable outcome in controller 
design. Moreover, one might argue that whichever of the above approxi- 
mations we take, we are in effect generating step like controls. Hence 
all our considerations which led us to proposition 2.5 are essentially _ 
violated: we are still faced with a problem of rate saturation and 

may even have worsened the situation by introducing approximations. 

On closer examination, however, all these problems are non- 
existent, if we do exactly what proposition 2.5 says - its actual goal 
is a continuous time coirpensator which generates a sequence of ramp functions 
as the control signals. Specifically, 
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x(t) = Ax(t) + B u(t) (2.108) 

u(1:) = u(kT) + (t-kT)u(kT) (2.109) 

for kT < t < (k+l)T. Note that ii(kT) , u(kT) only depend on the values 
of u(0) , x{0) , x(T),..., x(kT) , and j^(0) , # . . . (k-1) T) 

(equations (2.98) , (2.95)) . These computations can all be carried out by 
a digital computer in a sampled data fashion, provided one has a device 
to generate u(t) , given in equation (2.109) (e.g. an integrator) . Using 
equations (2.95), (2.73), (2.74) and (2.80), this can be written as 

u(t) = u(kT) + (t - kT) [K x(kT) + ^u(kT) + (K - La"^B) (C af } 

X A X \J 

( 2 . 110 ) 

There does not exist any non-causality problem at all. 

2.5 Final Remarks 

In this chapter We have presented several deterministic designs 
of compensators that can be used to control a linear time -invariant 
system to track a step input in continuous time, and discrete time. 

In the next two chapters, we shall concern ourselves with the formula- 
tion of the linearized longitudinal dynamics of the F-8c aircraft and 
the C* handling qualities criterion, and then formulate a control 
problem using some of the results obtained in this chapter.^ 

Comparisons among these various designs are made in Chapter 5, 
via simulations. It is shown there that the approximate schemes pre- 
sented in Section 2.3 are in general less effective than the scheme in 
Section 2.2. 

Before we close this chapter, let us compare briefly the servo 
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controller designs by T.A.S.C. [4] with these results. 

Firstly let it be said that in places where there are similarities,. 
T.A.S.C. 's approach looks more similar to [2], while ours is closer to [1] . 
In particular, instead of penalizing ^ directly, they penalize 

3 ^ - x(t) , which is clearly equivalent as shown in [1] or [2] or even in 
the proof of proposition 2.2. Also [2] shows that the techniques can be 
extended to the case when one puts an additional penalty term on 
u(t) - 1 ^, which is often done in [4] -. Furthermore, if one cares to 
do the algebra, one can easily see that the invertibility 

by T.A.S.C. (and [2]) is necessary for 
the assumptions (i) and (iii) in proposition 2.1, and the same holds 
for the corresponding results, too. Having said this, we shall consider 
the above to be equivalent conditions. 

The basic similarities between the two approaches are the issues of 
P-I ccntroller designs. Indeed, the derivations and results of T.A.S.C. 's 
continuous time and discrete time designs are the same as our results for 
system I and II respectively , except for the choice of u(0). They consider 
also issues of E-I-I controller, which we do not touch at all. 

The discretization scheme of the continuous time system and the 
cost penalty presented in [4] is completely the same as that derived in 
[7], which we referenced. Hence T.A.S.C. considered the case that u(t) 
is constrained to be step functions, and no detailed consideration has 
been given to problems of rate saturation. We constrain u(t) to be ramp 
functions, and derive the corresponding optimal control law, and discuss 
also the apparent non-causality problem. 


assumption of the matrix “ 


'A B' 

LX 

Ls. 2j 


CHAPTER III 

F-8C AIRCRAFT LONGITUDINAL DYNAMICS 
AND THE C* CRITERION 


i In this chapter, we describe briefly a linearized deterministic 
model and a C* handling qualities criterion of the F-8C aircraft longitu- 
dinal dynamics. Then we formulate a control problem (of the tracking type) 
on this linear model, using the C* criterion, in such a way that we can 
apply the results developed in Chapter II. 


3.1 Linear Equations of the F-8C Aircraft Longitudinal Dynamics 


The F-8C aircraft is a nonlinear system whose dynamics depends on a 
set of parameters /including its Mach number, altitude and dynamic pres- 
sure. Each such set of parameter defines an equilibrium flight condition 


forj the aircraft. J. Gera [141 linearized the nonlinear longitudinal 


dynamics about various selected flight conditions to arrive iat linear 
equations for the longitudinal motion of the F-8C aircraft. 

In this report. We elect to use sixteen flight conditions , namely, 
fli^t conditions 5-20, from ! [14], which represent a clean wing-down con- 
figuration of the vehicle. Note that this is the same set of data used 
by [8] . Indeed one of the objectives of this research is to investigate 
the feeusibility of incorporating the resulting controller design of this 
report to the MMAC scheme [8]. 
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In the linear model, there are four state variables inherent from 
the longitudinal dynamics of the aircraft: 


State Variables 

Units 

q(t) 

pitch rate 

rad/sec 

v(t) 

velocity error 

ft/sec 

a(t) 

angle of attack (from 
trimmed value) 

rad. 

e(t) 

pitch attitude (from 
trimmed value) 

rad. 


and 6 (t) is the perturbed elevator angle (in rad.), from its trimmed 
e 

value, driving these states. The linearized dynamics, for each flight 
condition, are given by 









q(t) 

-r; 

^11 • 


q(t) 

d 

v(t) 

— 



V(t) 

dt 

a(t) 


a * 
31 

^34 

a(t) 


0(t) 


1 0 0 0 1 

6(t) 








^NG 


5 (t) 
e 


(3.1) 


and in the servo-system, the actuator dynamics are represented by 


6g(t) = -12 6^(t) + 12 6g^(t) (3.2) 

where is the commanded elevator angle (in rad.). Furthermore, 

we include in the states a wind disturbance state w(t) (in rad.) . Choosing 

<5 (t) as a control, the conplete dynamics can now be written as, 

. ec - 

x(t) = A x(t) + B u(t) (3.3) 
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where 


x(t) = 


q(t) 

V(t) 

a(t) 

e(t) 

(S (t) 
e 

w (t) 


(3.4) 


u(t) = 

ec 


(3.5) 


A, B are 6x6, 6x1 matrices respectively, given by 


A = 


B = 


• 

A 

B 

^13 ' 

-LNG 

-LNG 

®23 

(4x4) 

(4x1) 

^33 

j'~' 


0 

0 

“12 

0 

(2x4) 

L 

■ 0 

0) 


0 

(4x1) 

I 

! 12 ;-- 

0 . 


(3.6) 


(3.7) 


where oo is the windpole depending on the flight condition, given by 


tji - ”2 

P 


{3.8) 


where V is the velocity of the aircraft, i.e. 
o 


V = Mach no. X velocity of sound (alt.) 


(3.9) 
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f , 

! 2500 ft. 
200 ft. 


if alt. > 2500 ft. 
if alt. = 0 


(3.10) 


We assijme that our sensor takes two tneasurements , one on pitch 
rate, q(t) , and the other on normal acceleration, a.^{t) , which is given, 
about equilibrium, (positive in the upward direction) by 

V ■ 

a^(t) (q(t) - a(t)) V (3.11) 

z g 

-2 

where g is the acceleration due to gravity (32.2 ft, sec ). 

Equations (3.11) and (3.1), and the fact that a^^^ = 1, = 0 for all 

the flight conditions , give us a linear combination of the state variables, 

i.e. , 


, , V : ■ 

a (t) = -^ (-a,^ V(t) - a _ a(t) - b 6^(t)) (3.12) 

z 9 ® 


Therefore, in the form of an observation matrix equation, we 
write ■: 



y(t) 

= C x(t) 



where C is a 2 

X 6; 

matrix given by 




.1 

0 0 

0 

0 0 


0 

V32 V 33 

g g 

0 

V h- 

- ° 3 Q 

9 


and 


(3.13) 


(3.14) 
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x(t) = 


q(t) 




(3.15) 


This completes a brief deterministic description of the linearized 
P-8G aircraft longitudinal dynamics. 


3.2 F-8C Aircraft C* Criterion 

In this section we summarize and discuss briefly the significance 
of a C* handling qualities criterion in the longitudinal dynamics of a F-8C 
aircraft. 

In today's technology, the implementation of flight control systems 
in an aircraft has been widespread. In order for an engineer to make 
sound judgement on the merits or short-comings of a control design, he 
needs a well defined criterion on the handling qualities of the aircraft. 

A proper choice of such criterion is clearly essential to demonstrate the 
effectiveness of a good flight control system. 

It used to be a. common practice in classical designs to use a "short 

period damping ratio" type of criterion. This is becoming out-of-date for 


the follcwing , reasdns j 

(I) it assimes that only one variable is predomiriantly sensed 

! 

by the pilot, namely, ithe normal acceleration; 

C2) it assumes that the short period response of tiie system may 
be represented by a linear second order system; 

(3) it does not give an exact picture of what the aircraft is 


really doing in the t reins ient period. 

For a high performance aircraft such as the F-8C in which maneuvers 
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are frequent due to its agility, this is certainly not -satisfactory . It 
has been observed that normal acceleration, as well as pitch attitude and 
pitch rate, are all important cues to a pilot, and their relative impor- 
tance depends on the velocity of the aircraft [10] . Specifically, normal 
acceleration is more important when the velocity is high, and the pitch 
variables predominate when it is low. Moreover, an aircraft is more com- 
plex than any simple representation by a second order linear system. 

Also it seems reasonable to rely on a criterion in the time domain, 
rather than one in the frequency domain, partly because of its ease to 
perceive, and partly because of the past success of modern control theory 

Thus it is intuitive that in order to develop a more satisfactory 
criterion of aircraft handling, we may want to first establish a time 
history envelope for one or more chosen paraineters of the aircraft dynam- 
ics. Such an idea is well illustrated in [10] , [11], by a so-called c* 
handling quality criterion, and indeed this is the step that we take. 

We *;^ive a brief formulation of this criterion in terms of the F-8C 
aircraft longitudinal dynamics below. 

A variable C*(t) is defined to be a linear combination of the normal 
acceleration, a^(t), and the pitch rate, q(t) , which is given by the 
following relation (in g’s) 

C*(t) = k. a (t) + k qAt) (3.16) 

■ ■■ I M Z ^ 

where k^ , k^ are some arbitrary pair of constants, and can be viewed as , 
the relative weightings of the importance of these two state variables 
in the pilot's sensing. Arbitrarily, k^ is chosen to be unity# and 
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correspondingly , is given the value 

, cross-over velocity (ft/sec) 

^2 “ 2 {3.17; 

g (ft/sec ) 

where the cross-over velocity is the aircraft speed at which the normal 
acceleration and the pitch rate give equal cues to the pilot. Approximately 
is thus found to be equal to 10 with respect to the F-8C aircraft 
dynamics, or equivalently, 

C*(t) = a (t) + 10 q(t) (in g's) (3.18) 

z , . 

Figure 3.1 shows a time history envelope for the C* parameters thus 
defined, representing an acceptable range of the C* transient response to 
a step input from the pilot. 

3.3 Continuous Time Formulation of a C* Tracking Design 

Thus far we have given almost all the preliminaries required for a 
deterministic compensator design, such that the C* transient response of 
the P-8G aircraft to a step input is made to fall within the C* time 
history envelope, e.g. along the middle dotted line of the envelope in 
Fig. 3.1. In what follows, our task is to give a physical motivation 
for the overall design, and then discuss the varioxis theoretical issues. 

In the next section we shall extrapolate these ideas into the discrete 
time designs for the F-8C 'aircraft. »■ 

When a pilot attempts to control an aircraft, most probably he applies 
an inpluse like action to the command stick, inducing a step like shift 
in the stick position. This in turn is transformed into a pilot's desired 
input value to the dyneunics by some servo-mechanisms. With regard to the 
C* criterion, we imagine that such a step input thus generated would be a 
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O 



TIME ( SECONDS) 


Fig. 3.1 C* Time His^ory Envelope (Regions 1 and IV ore occepfable for the C* response) 
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desired steady state C* response that the rilot wants to have the 
short period response of the aircraft to frllow. We summarize these in 
Fig. 3.2. 

Note that in this model we have assuned that a G* desired value is 
generated by the pilot at time t^, represented by , and there 
is no time delay in the pilot's reflexes aid motion, nor in the servo- 
mechanisms of the aircraft. Prom now on, ve shall assume implicitly 
this part of the model, and only consider the desired C* response, as 
a step input to the aircraft at time t^. vithout loss of generality, we 

take t = 6. 
o 

• A 

The formulation of the control problsn that we are aiming at is now 
clear. Specifically, we are given the lirear longitudinal dynamics of the 
F-8C aircraft, described by the equations 3.3) - (3.11) in section 3.1, 
and we want to find an optimal control function 

u(t) = 6 (t) t e [0, w) (3.19) 

ec 

such that the following performance index is minimized: i 
a(u) =/ [(CJ 

where Q, R are some chosen non-negative nunbers. 

We divide the specification of this problem into two parts: 

(i) Statement 

The deterministic descriptiois of the dynamical system 


C*(t)fQ + (t)R I dt 
ec 


■] 


(3.20) 
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Fig. 3.2 C* Pilof-AircToff InteracHon Model 
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x(t) = A x(t) + B u(t) 


C*(t) = C* x(t) 


(3.22) 


where A,B, x(t ) , u(t) are described by equations (3.4) , (3.5) and (3 .7) , (3.8) 


In particular, is a 1 x 6 matrix, given by 


C* = 


V 

T ^32 


V 

“ a 
g 33 


• -f‘. •] 


(3.23) 


In addition, we are given a desired steady state value of C*(t) , namely 

C*, which serves as an input (step at t=0) to the system. We want to 
a 

find an optimal, everywhere continuous control u(t) such that this per- 


formamce index is minimized. 




(C* - c*(t)) Q + 


where Q, R are positive real numbers. 


(ii) Solution 


•2 

U (t)R 


(3.24) 


It is easy to see that this belongs to the class of problems 
examined in proposition 2.1, with scalar output and control, if we assume 


i) A is non-singular; 

ii) [A,B, C^] is minimal; 

iii) C*A is invertible. 

For the purpose of this chapter, let us simply assime them to hold and 
we shall see in Chapter V that all the flight conditions satisfy ii) and 


iii) , and except for a few of them, i) is also satisfied. The exceptions 
are due to mmerical approximations in the data. 

By proposition 2.1, we are then able to find a scalar constant L and 
a 1 X 6 matrix N such that 

u(t) = lJ (C*-C*(t)) dT + N{x(t) - x(0)) + u(0) (3.25) 

0 

is the optimal control solution. Fig. 3.3 shows the structure of the con- 

n 

CT V (X fc) s w 

trolled aircraft, where N,N,N,N,N ,N are entries in the N matrix. 

3.4 Discrete Time Formulation of the C* Tracking Design 

Recall that in Chapter II, we discuss two schemes of discretizing 
the continuous time dynamics , i.e. system II and III. Here, for each of 
them, we outline a formulation of a control problem, on the F-8C aircraft 
longitudinal dynamics, which is equivalent to the one just given for con- 
tinuous time. Because of the possibility of too many rediandant discussions 
with Chapter II, and the many similarities to the continuous time problem, 
much of the detail here is omitted. For simplicity of discussions, we 
assume conditions such as controllability, observability and invertibility 
of matrices implicitly whenever they are needed. 

1) Discretization using Scheme (1) 

The continuous time linearized F~8C aircraft longitudinal dynamics, 
described in section 3.1, Ccih be discretized by constraining u(t) (i.e., 
6ec(tj) to be a constant over each sampling period T, i.e. for each k and 
t e [kT, (k+l)T), 
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u(t) = u(kT) = u 


(3.26) 


with the application of proposition 2.3. Specifically, we get a deter- 
ministic discrete time realization 


where 


“ic+1 -^-Tc -d k 


G* = C*^ 
k k 


(3.27) 

(3.28) 


A at 
-d 


2d 


= / 


dT B 


(3.29) 

(3.30) 


Suppose in addition, we are given a steady state desired value of C*(t) , 
say C^, and we are required to minimize the performance index 


J(u) - £ KCJ - 4- (u 

js=0 ■' 


(3.31) 


where I Q^, are positive real numbers. Then, by applying proposition 2.2 

we can find a scalar L,, and a 1 x 6 matrix N such that the optimal dis- 

d “d 

Crete time control law is given by 

■ ■ k-1' 

% = L , 2 (C* - C*) + N,(x. - X ) + u 
k d d j -d He -o o 


(3.32) 


Instead of considering a performance index such as equation (3.32) , 
one might want to use one of the two approximations of the continuous time 
desvign, discussed in Section 2.3. Then the control law is given by 
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k-1 

j=0 


(3.33) 


where L„, is a scalar, N_ is a 1x6 matrix computed by using either 

u D 

equations (2.56) and (2.57) or equations (2.63) and (2.64)., 

2) Discretization using Scheme (2) 

Alternatively, we may discretize the continuous time problem by 
keeping <Sg^,(t) to be a constant over each sampling period T, i.e. for 
each k and t e [kT, (k+l)T), 


u(t) = u(kT) = Uj^ 


(3.34) 


Then by proposition 2.4, we obtain the following discrete time realization 
for the deterministic c* dynamics. 


*x * 

A +1 


X' 




3 t+i. 


u 

L k J 


+ D 




(3.35) 


where 


C* = [C* 0] 


\ 


(3.36) 




^ 2d 


(3.37) 
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D' = 


T S , 

f J e-^ dx ds B 
o o 


(3.38) 


Matrices ^ are given by equations (3.38), (3.39). Moreover, we arrive 
at a performance index equivalent to that for continuous time (equation 
(3.25)) , 

J = Z 2.% - 


where 


2 


M 


= C*"C* e-^ 


dt 

Jt ^s 


= Q 7 e- ^ C*"C* if f e-^ 


dT ds]B dt 


(3.40) 


(3.41) 


o o 


and 


R = PT 


X, = a“^b (c*a”^b) “^C* 
“d — d 


(3.42) 


(3.43) 


Applying proposition 2.5, the optimal sequence such that this perfor- 
mance is minimized is given by 

'■k-l . 

u = £ y] (c* - C*) + ?i[(jc - X ) + u 
k ^ d 3 -He -D o 


(3.44) 


where L is a; scalar and N is a 1 x 6 matrix found by the procedures given 
in that proposition. Along the line, we also obtain a lx 6 matrix and 
a scalar such that the optimal control rate for t e [kT, (k+l)T) is 
given, 'by 
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i(t) ■■= u(kT) = ^ = K^(Xj^ - + l2^"k ^ V 

where 

u = -(C*a”^B)”^C* 
d ; ;d 


(3.45) 


(3.46) 


and the control generated should be a sequence of continuous time ramp func 
tions, to be applied to the continuous time system, given by 


u(t) = Uj^ + (t - kT)Uj^ (3.47) 

Thus we have completed the discussions on full state control tracking 
designs for the F-8C aircraft longitudinal dynamics, using the C* handling 
qualities criterion. In Chapter IV, we present a subpptimal design, by 
reducing the number of states, which we compare with in Chapter V on the 
basis of numerical results. However, we only consider system I (continuous 
time) and system II (discrete time) , and leave system III for future 


research. 


CHAPTER IV 


A SUBOPTIMAL REDUCED STATE DESIGN OF THE 
E-8C AIRCRAFT LONGITUDINAL DYNAMICS 


4,1 Preliminairy Remarks 

. < -I ■ ' i . 

In [8], it is sketchily pointed out that a reduced state design can 
be an! extremely good approximation to the transient behavior of the full 
statej F-8C aircraft longitudinal dynamics, particularly in the stiboptimal 
regulator design. The purpose of this chapter is to exploit this idea 
in the C* design formulated in Chapter III. Treatment here includes 
a brief account of the desirability of such a reduced state design, and 
some detail of the design equations. We comp<ire and contrast this design 
with the optimal one along the ^Aeoretical issues. Numerical comparisons 
are given in the next chapter. 

Here we derive how one can go from a full state model to a reduced 
state model in continuous time by making approximations and linear trans- 
formations J To obtain a reduced state model in discrete time, we remark 
that this can be accomplished by discfetizing the continuous time reduced 


state model using either one of the methods given in Chapter II, as 


we have done For the full state model in Chapter ml Bearing this in 

, ! >■ i f ^ ^ ■ 

mxnd, we restrict ourselves to discussions in the continuous time case 
only for the rest of this chapter, since the discussions in Chapter III 


on discrete time systems serve the purpose equally well here, provided that 
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one remembers to deal with a reduced state system instead of an optimal 
one . 

Note that notations used in Chapter ill are also used here with 
the same meanings? no attempt is made to redefine them. 

4.2 Reasons for Studying the Reduced State Model 

Prom a study of the components of the system matrices. A, of the 
different flight conditions, it is noted that two of the states, ncunely, 
the velocity error, V, and the pitch attitude, 6, seem to have very smell 
effect on the short term transient responses of the F-8C aircraft. Thisi 
is further confirmed by an analysis of the eigenvalues and eigenvectors 
of the system matrices. Apparently, we can eliminate these two variables 
from the model design, without degrading much of the short period perfor- 
mance of the aircraft. 

Such an intuitive idea is pursued in [8], in the design of a siabop- 
timal feedback regulator. It is noted that the feedback gains on V and 
e are very small in the full state regulator design. Indeed, feedback 

gains on the other state variables, i.e. q, a, 6 , 6 , w, are approximately 

: e , ' ec ' ■ ■ 

the same as those from a reduced state model design, which uses only these 
five variables as states. Different plots of closed-loop poles, and gains, 
plus a very detailed analysis of simulations, show that the transient 
responses of the aircraft in both cases appear to be very similar. 

It is intuitive that if a reduced state design is acceptable, we may 
expect that the overall engineering design (implementations, etc.) can be 
simplified. Eliminating the state varicibles V and 0 from our 
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full state model in Chapter III, we 
dynamical representation; i 


_d_ 

dt 



q^t) 



a(t) : 


<s (t) 



e 


L 

w(t) 

. 



31 

0 

0 


^13 


^33 


0 

-12 

0 

0 


arrive at the following (deterministic) 




i 

r "1 



- 


q(t) 



0 


a(t) 



0 




- 



6 (t) 


12 


e 




w(t) 


0 

. . 


6 (t) 

ec 


(4.1) 


We take a step further from this by replacing the state variable ot(t) 
by a new (but familiar) state variable, a (t) , the normal acceleration. 

Z ' 

This can be easily accomplished by techniques such as matrix transforma- 
tion, giving us 


q(t) 


^11 ^ 

^o ^33 

, 13 

*13 


q(t) 

a^(t) 

0 ■ 


"®33 g 

^33 

3 g 

g 33 


a^(t) 

«^(t) 


0 

0 

-12 ■ 

0 



w(t) 


0 

0 

0 

0) 


w(t) 


'i ■■ 'i 0 1 

*^o 
-12 bL -^ 

1 . . ,3 g 
1....12 , 

0 




To simplify notations, we write this as 


^(t) = ^ x^(t) + u(t) 


(4.2) 


(4.3) 
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Hence the C* parameter is given by 

' C*(t) = [10 1 0 0]x (t) (4.4) 



This tremsformation is justified for the following reasons: 

(i) Computations, whether on-line or off-line, are reduced; 

(ii) Since ve make measurements on the pitch rate and normal accel- 
eration, and the C* parameter depends only on these two states, therefore 
such a transformation provides us the convenience of having a constant 
observation matrix C^, cind a constant matrix £* for all flight conditions; 

(iii) We can form C*(t) directly from the measurements. 

We reinark further that except for the use of equations (4. 3) and 
(4.4) , a C* tracking design on this model can use the same performcince 

index used in the full state system design, i.e. 

-00 

J = I [(C* - C*(t))^ Q + u^(t)R] dt (4.5) 

*^o 

Using this, and applying proposition 2.1, we thus obtain a scalar L^, and 
and a 1 x 4 constant matrix 

N = [n n n^ n 1 (4.6) 

z e 

such that a sxiboptimal deterministic design is achieved r with a control 
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law given by 


uCt) ^ f 


(C* - C*(T))dT + Nj.(Xj.(t) - ^( 0 )) + u(0) (4.7) 


Pig. 4.1 shows the structure of this reduced state controller, and 

a 

is to be compared with Fig. 3. -3 for the full state design, where ' 

w . . 

N_ , N are entries in the N matrix. Equivalently, we can write this as 


u(t) = L 


r 


(C* - C^(T))dT + N_<x(t) - X(0)) + u(0) 

d I “ 


(4.8) 


where 


(4.9) 


and x(*) is the full state vector. 

We shall analyze the effect of such a stiboptimal design on the 
full state aircraft model, by studying some numerical examples in 
Chapter V. 


Nf = Nr 


(2 X 6) 


(2x4); (2 X 2) 










CHAPTER V 


NUMERICAL EXAMPLES 

5.1 Overview 

In this chapter, we present some nvunerical results on the F-8C air- 
craft linearized longitudinal dynamics that illustrate our work on the 
deterministic design, both in continuous time and in discrete time. We 
selected to use sixteen of the flight conditions (5-20) from Gera [14], 
which represent linearized dynamics of the F-8C aircraft in clean wing dow‘.i 
configurations. Table 5-1 shows each of the sixteen flight conditions 
(altitude, Mach no., dynamic pressure). Throughout this chapter, we em- 
phasize the fact that the full state model of the F-8C aircraft (see Chapter 
III) is very well approximated by the reduced state model discussed in 
Chapter IV. Thus a major part of this chapter contains results for the 
reduced state designs, although brief notes of comparison with the full 
state model are given along the discussions. 

We remark that in linear simulations on a digital computer, one in 
general does not encounter the problem of rate saturation of the control 
vfhiv'n is discussed in some detail in the later part of Chapter II. There- 
fore, for the purposes of this report, we only use System II (See proposi- 
tions 2.2, 2.3) and leave aside System III (See proposition 2.4). Also 
in Section 2.3, we presented two discrete time designs which make use of 
and represent (approximate) discretizations of the continuous time design. 
We shall not give any numerical results for them, but instead we point out 
here that simulated results have shown that both of them are less effective 

■' - 66 - ■ 
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Flight Condition Altitude Mach. No. 


Dynamic Pressure (PSF) 


5 

sea 

level 

0.3 

133.2 

6 

If 

II 

0.53 

416 

7 

II 

II 

0.7 

726 

8 

It 


0.86 

1098 

9 

II 


1.0 

1480 

10 

20 , 

000 ft. 

0.4 

109 

11 

II 


0.6 

245 

12 

II 

II 

0.8 

434 

13 

II 

It 

0.9 

550 

14 

II 


1.2 

978 

15 

40, 

000 ft. 

0.7 

135 

16 

II 


0.8 

176 

17 

II 


1.9 

223 

18 

II 

II 

1.2 

397 

19 

II 


1.4 

537 

20 

II 


1.6 

703 


Table 5.1 Flight Condition Parameters 
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than the discrete time design derived directly from the dynamical system 
(See propositions 2.2, 2.3), though the first order approximation method 
gives rather satisfactory results. Indeed the one obtained by pole-alloca- 
tion is far from satisfactory and should not be used in practice. 

There are five different designs to be investigated in this chapter, 
including two continuous time designs, and three discrete time designs. 
However, only the ones that deal with the reduced state model will be dis- 
cussed in detail. We summarize them in Tcible 5.2a. 

We selected three’ of the sixteen flight conditions for simulations in 
order to demonstrate the effects of mismatching in the C*-controlled flight 
trajectory. Namely, they are flight conditions 7, 11 and 12, representing 
variations in dynamic pressures, with flight condition 7 being the highest, 
£Uid 11, the lowest. We assvime that the aircraft is in flight condition 12,. 
while control gains from one of these three flight conditions are used. 

We remark that all of these flight conditions are subsonic, where the 
C*-theory plays the most significant role. We shall assume that initially 
all the states are zero and the desired C* steady state value is a step of 
1 g. Also u(0) = 6ec (0) is chosen to be zero. We summarize them in 
Table 5.2b. 

From now on, for simplicity purposes, with reference to Tables 5.2 
a and b, we shall refer to the simulation i, l<i£3, for the design 
n = A,B,C,D,E, a^ the simulation J2-i. We remark that simulations 1 and 2 
represent mismatches in gain scheduling. Where appropriate, we shall attach 
plots of simulations to the discussions. All the simulations are for 6.25 



Gains 


Design 

A 

B 

C 

D 

E 


Time Set 


continuous 


discrete 


Simulated Model 

reduced state 
full state 
reduced state 
full state 
full state 


from reduced state model 


II tl 


II It 


from full state model 


Simulations 

A 

B 

C 

D 

E 


Table 5.2a Summar^^ of the Designs Simulated 


Simulation 

Flight 

Condition 

Control 
Gains (F/C) 

Initial Conditions 
of States 

C*d(g's) 

1 

12 

7 

i all zero 

1 

2 

12 

11 

(i) < " " 

1 

3 

12 

12 

1 II II 

1 


Table 5.2b Simmary of the Simulations for Each Design in Table 5.2a 
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seconds r with a time step of 0.125 second, which is the sampling period 
for discretization. 

For each of the two reduced state designs (designs A and C) , we show 
a iriismatch stability table of the sixteen flight conditions, which indicates 
whether the flight woyld be stable when control gains computed for flight 
condition i (5‘^i£20) , are scheduled to be used in controlling the aircraft 
in flight condition j (5^^20). These tables, together with the simulations, 
form the basis for the mismatch stability analysis in this report. 

We arrange the rest of the chapter as follows: 

Section 5.2 deals with the continuous time designs (designs A and B) . 

There brief notes are given for the choice of the weightings in the cost 
function, and mismatching of the flight conditions. A very detailed analysis 
for the reduced state design is presented. Control gains are computed for 
the reduced state and the full state models. The same is done in Section 
5.3 ifor the discrete time designs (designs C and D) . In addition, a full 
state model with full state control gains (design E) is also discussed. 

In Section 5.4, we draw some conclusions and add some remarks. 

For simplicity and reference purposes, all notations used in this chap- 
ter are the same as those in Chapters III and IV, unless otherwise stated. 

It should be noted that use has been made of Sandell and Athens [16] 
and Kleinman [17], in the computing work done for this report. 

5.2 Continuous Time Designs 

The first step that one should take is to choose the appropriate 
weightings; Q and R, for the cost function given by equation (3.25). 
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In doing so, we pick R from a consideration of the rate saturation of the 
actuator dynamics; 

1 (5.1) 


R = 


ec 


max 


where 


6 = 0.435 rad/sec (25 deg/sec) 

max 


(5.2) 


Then we chose Q in such a way that the short period C* responses of the 
closed loop dynamics fall more or less in the middle of the G* time his- 
tory envelope shown in Fig. 3.1, i.e., along the dotted line. We found 


out, by trial and error, that a good choice for Q is 
Q = 1 (g"^) 


Figure 5.1 shows representative C* responses. 


(5.3) 


(i) Design A' 

Using this pair of values of Q, eq, (5.3)> and R, eqs. (5,1) and (5.2) 
one can then proceed with the (continuous time) design discussed in Chapter 


IV for the reduced state 


model. That is given equations (4.3), (4.4), 


and (4.5), one is to find the control 5 (t) given by equation (4.7). 

It is observed that, for all flight conditions , the matrices , 

B^ > C * satisfy all the three assumptions required by proposition 2.1. We 

■ I. - : 

can therefore apply the design procedures given by this proposition and ob- 
tain the gains L^, NT directly. 

Figure 5.2 shows a plot of the resultant dominant coitplex poles of 
the closed loop system for all the flight conditions. The magnitude of 



NORMALIZED 



Fig. 5J Illustrations ofhow Choice Of Q is Made. Shown Above are 
Simulations For 3 Values of Q on Flight Condition 13. In 
General, Larger Values of Q Induce Faster Modes of C* 
Responses. Q = 1 is Chosen for this Report 
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the complex pole increases as the dyneunic pressure increases, at each alti- 
tude. These, of course, correspond to the short period dynamics. For 
the purposes of comparing this C* design with some other traditional handling 
criteria of aircrafts, we draw in three straight lines to show the range 
of dancing ratios through these flight conditions. One such criterion 
is the Cornell Aeronautical Laboratory "thvunbprint" [10]. This "thumbprint" 
represents combinations of short period natural frequency and damping 

ratio (^) that were preferred by pilots of several variable stability 
aircrafts, and is valid on the assumption that the predominant variable 

sensed by the pilot is normal acceleration, a . Figure 5.3 shows such a 

. ' z 

thumbprint plot together with three vertical lines corresponding to the 
damping ratios in Fig. 5.2. Comparing these two figures, we see that all 
the flight conditions lie more or less within the thumbprint, except flight 
conditions 9 and 14 which represent transition behavior between subsonic 
and supersonic flights, Tlie fact that flight conditions 18-20 lie close to 
the boundary is not surprising since they are supersonip. Most importcint 
of all, we observe that flight conditions 5, 10 and 15 lie very close to 
the boundary. This is explained by tiie fact that they have low Mach numbers, 
emd indeed, the pitch variables are more important cues to the pilot. 

Tcib'le 5.3 is a mismatch stability table, in the sense that the 
aircraft is in flight condition i, 5^1^20, but one is using control gains 
calculated from flight condition j, 5<j<2o, Specifically, the closed loop 
dynamics are 

x(t) * A x(t) + B u(t) (5.4) 

” ■ ,“T ■ ■ ■ ■ —IT 

. ■ ^ 1 i. 




STiiEILITY SUMMARY TABLE 


CONTROLLER 


TRUE 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 . 

16 

17 

18 

19 

20 

FC 

















5 

♦ 


♦ . 

♦ 

4c 

4c 

4c-. 

4c 

4c 

4c - 

4= 

4c 

4c 

• 4c 

„ ,4c^ 

4c 

6 

♦ 

♦ 


4c 

4c 

4t . 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4= 

4c 

4t 

7 

♦ 


' * 

♦ 

4c 

4c 

4c 

4c 

. * ' 

4c 

4c 

4c 

-4c 

'4c 

; 4c 

4c 

8 ' 

♦ 


4c 

♦ 

♦ 

4c 

4c 


4t' 

4c 

4c 

-4c 

4c 

4c 

4c 

4c 

■ ;9 - 

♦. . 

♦ ■ 

' 4> 

4c 

4c 

4t 

4c 

4c 

4c 

4c 

4c 

4c 

4c 


4c 

4: 

ro 


♦ 

♦ 

♦ 

4c 

4c 

■ ♦ 

4c 

4c 

4c 

: 4c 

4c 

4c 

4c 

4c 

4c 

11 

■ 4c 


, ♦ 


4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 


■4c 

4c 

12 

♦ 


♦ 

4! 

* 

4c 

4c 

4c 

4c 

4c 

^ * 

4c 

4c 

4c 

4c 

4c 

13 

♦ 


♦ 

4c 

♦ 

■4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

1 4 

■ '♦ ; 


4« 

4c 

4c 

4c ' 

4c 

4c 

4c 

4c 

4C: , 

4c 

4c- 

4c 

4c 

4c 

15 

■ ♦ 

♦ 

♦ 

. ^ 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

: 4c 

4c 

4c 

4c 

4c 

16 

* 

♦ 

♦ 

4c 

4c 

4c 

4c 

4c 

4c 

* 

4c' ; 

' 4c 

4c 

4c 

4c 

4= 

17 

* 


♦ 

♦ 

4< 

4< 

4c 

4c 

4c 

4c 

‘ 4c 

4c 

♦ 

4c 

4c 

4c 

18 


♦ 


4c 

4c 

4c 

4c 

4c 

4c 

4t 

- 4c -■ 

4c 

4c 

4c . 

♦ 

4= 

19 

♦ 

♦ 

* 

4c 

4c 

4c 

4c ■ 

4c 

4c 

4c 

4c 

4c 

4t 

4c 

4c 

4c 

2C 

* 

' ■* 

♦ 

* 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4c 

4: 


n=[JNSTABLE 
*= STABLE 


Tcdsle 5.3 Mismatch Stability Table for Design A 

(All mismatches are stable.) 


-77- 


u(t) = L (C* - C*(x))dT + N (x(t) - x(0)) + u(0) (5.5) 

r . I d ~~x . — ’ 

where A . B are A , B matrices o£ flight condition i, and L and 
— r. — r. ~r -X • — r. 

11 3 

N are control gains from flight condition j . 

“X .... 

3 

It is seen from Table 5.3 that all mismatches give stable flight 

for the system. Also we can see from the simulations A-1 to A- 3 (Figures 

5.4 - 5.9) that the C* and 6 responses are in general very satisfactory f 

e 

c 

with acceptable control rates. When there is not any mismatch (i.e., i,j=12), 
the system gives the best response, as is expected. .Observe how nicely 
the C* respoKise fits into the C* envelope in simulation A-2 (Figures 5-6, 5-7) . 
Mismatches give a little slower response, and more overshoot. 

We attach also plots of the L and N gains against Mach numbers and 
dyncunic pressures (psf) (Figures 5-10 - 5-19) . Notice that in each plot 
the subsonic flight conditions are readily divided into three groups according 
to altitude, and the supersonic flight conditions form a single group by 
themselves. 

(ii) Design B 

In Chapter IV it was discussed how one can apply the reduced state 
control gains computed in design A to a full state model, in equations 
(4,8) and (4.9). We study such a design strategy by doing simulations. 

It is seen that the responses in all cases are almost identical to 
the corresponding ones for design A. Also we remark that this design gives 
stable flights for all mismatches. This suggests very strongly that the 
reduced state design is a good design for the system, at least as far as 
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Simulation A-3; C*-response to a unit step C*^ input; match (12-12) 
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continuous time is concerned. As an example, we attach Figures 5.20, 5.21 
which show simulation B-3 of the matched flight condition 12. 


5,3 Discrete Time Designs 

To pick appropriate values for the weightings Q^, R^, in the cost 
function given by equation (3.39), one might want to do it in such a 
way that this cost fiinction approximates that for continuous time, 
i.e., equation (3.21). We elect to use the following values: 


= tq 

= R/T 


(5.6) 

(5.7) 

where T is the sampling period, which is 1/8 sec. in our design, and 
the choice of Q and R is given by equations (5.1)-(5.3). This was done 


because we can use the approximations 
u( (k+l)T)-u(kT) 


I 


u(kT)~ 
kT 


f(T)dX~ T f(jT) 

0 j=0 


(5.8) 

(5.9) 


for any real continuous time function f(*) and hence (C*^ -€*(•)) Q. 
Thus, we use the following weighting parameters in Eq. (3.39) 



1 

8 

8 X 0.435 


(5.10) 

(5.11) 


(i) Design C 

Using the above values of and R^, and the discretized reduced 
state model obtained by applying the discretization scheme in propo- 
sition 2.3 to the continuous time reduced state system given by equations 
(4.3) , (4.4) , we can follow the design procedures in proposition 2.3 




(rad. ) 
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to find control gains L , and ^ . That is, we may write, 


X ((k+l)T) = X (kT) + u(kT) 

~~X ~Q —X —U 

r r 

k-1 


(5.12) 


u(kT) = L (C* - C*(jT)) + tL (X (kT) - X (0))+ U(0) (5.13) 

d ,M d -nd -r -r 

r 2=0 r 


Indeed the assumptions required by proposition 2.2 are found to be 
satisfied by all the sixteen flight conditions. 

Considering mismatches, we form the dynamics 

X ((k+l)T) = ^ X (kT) + u(kT) (5.14) 

~~T —a — r —a 

r. r. 

X 1 


u(kT) = L- (C* - C* (It)) + ^ (X (kT) - X (0)) + u(0) 

d d —r — r 

r.x,=0 r. 

3 3 


(5.15) 


where 5 <x, 3< 20 and A, , B are A, ,B, matrxces of flxght condition i, 

— — “d — d —a — d 

r. r. r r 

11 

and L » Nj control gains L , N from flight condition j . 

Q d u — a 

r r 

A mismatch stability table for this design is shown in Table 5.4. 

Here instead of universal mismatch stability, as in the continuous 
time case (Table 5.3), there are about 10% of the mismatches that 
are unstable. In particular, most of them are mismatches of controllers 
from flight conditions that have low dynamic pressures to flight condi- 
tions of high Mach ntuiibers (and high dynamic pressures) . Of course 
the system is always stable whenever there is no mismatching. 

Looking at the simulations (Figures 5. 22-5. 27) , we observe that they 
all show very desirable responses, indeed they are almost the same as 
those for design A, although now we have a little more overshoot. The 
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Table 5.4 Mismatch Stability Table for Design C 

(e.g. if the aircraft is in FC 9 (row) and one use the controller of FC 7 
(column) r the closed loop aircraft will be stable; if one uses controller 
of FC 12 (column), the closed loop aircraft will he unstable.) 


C*-response (g's) 






-response (rad.) 



Fig. 5.23 Simulation C-1; <Sec-response to a C* input; Mismatch (12-7) 
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contrpls differ somewhat, with higher initial rates of change here, but 
the slopes are more or less the same. Again these mismatches seem to 
impose no problem at all in the overall control of the aircraft. 

We also attach plots of the L and N gains against Mach numbers 
and dynamic pressures (psf) (Figures 5.28-5.37). Notice that, similar 
to the continuous time case (design A) , in each plot, the subsonic 
flights are readily divided into three groups according to altitude, 
and the supersonic flight conditions form a single group by themselves. 


(ii) Design D 

Similar to the continuous time case, we can use the control gains 
from the reduced state design C, in the full state model. That is, we 
generate the control by the relation 

u(kT) = L (C*. - G*(jT)) ■+ N (x(kT) - x(0)) + u(0) (5.16) 

Y* fl ' ■- ’ r — — • 


'd j=0 


where 


^ = N 

■^d ^d 


2x6 

0 0 0 0 1 0 

0 0 0 0 0 1 


(5.17) 


The simulated results are almost identical to those in design C, 
and in fact, are very similar to the continuous time simulations (de- 
signs A and B) . 


(iii) Design E 

In doing a reduced state design, one is using less information 
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than what is available to him. One may expectv therefore, that a full 
state design is more accurate and perforins better than the reduced 
state design. The full state design is described by equations (3.35) 
to (3.40) . 

Con^aring the simulations of such a design, with design D, it is Idbserved that 
both the responses and controls almost coincide exactly in the two 
cases. As an exait^le. Figures 5.28-5.29 show simulation E-3, which 
is a matched simulation for flight condition 12. 

If one examines the gains, one ceui see that they are very close 
together for this design and the design D. 

5.4 Final Remarks 

In general, the reduced state design does seem to provide a 
very good approximation to the full state design. As far as responses 
are concerned, all designs (A-E) are very similar. This shows a strong 
correspondence between the continuous time designs and the discrete 
time designsL Higher initial rate of change of the control occurs 
for the discrete time designs, however. Nevertheless, they all seem 
to be acceptable. 

If one con^ares the plots of control gains for the continuous 
time reduced state design (design A) with those for the discrete time ^ 

reduced state design (design C), one can see that the N gains are in 

general rather close together, with the continuous time gains being 

somewhat smaller. Also the continuous time L gains are approximately 
one eight of those in discrete time. This accounts for much of the partial 
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success of the method of first order approximation of the continuous 
time design (although results in design C are better) . 

The above discussion holds, at least as far as the linearized 
system dynamics go. In particular, it confirms our theoretical results. 


sx 
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CHAPTER VI 

CONCLUSIONS AND SUGGESTIONS 

The approach taken in this report seems to provide satisfactory 

results on linear time-invariant dynamical systems, which are required 

• i 

to track step inputs, both from a theoretical and a practical point of 
view.: Theoretically, it is illustrated by the propositions given in 

Chapter II, where we give diigressions of both ;the continuous time and 
discrete time designs. Their applications to the Pr8C aircraft linear- 
ized llongitudinal dynamics:, ^ing a C* handling qualities criterion, 
provide usi a lot of physical’ insights, as we have observed from the 
control gains euid various simulations. 

Apparently the reduced state designs possess very high practical 

■ ' ; i 

values. Chapters III and IV illustrate in some detail how they can be 
implemented in the aircraft. It is believed that further study into 
this can be very fruitful. 

The two discretization schemes discussed in Chapter II, of which 

. i ' I - ■ ' ■ ■ ■ ■ 

one keeps” the control to be a constant over each sampling period, and the 

other, its rate, bear very inportant significance in physical interpre- 
tation, because of the possibility of rate saturation in the actuator 
dyneunics. We think the different controls computed by these two different 
schemes should generate different performances of the system. Since 
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scheme II constrains the control to be continuous, it is likely that it 
approximates the continuous time system better. Apparently from Chapi3r V, 
the continuous time design provides the best result. However, the dxrital 
requirements in implementation limit its usefulness. On the other haid, 
the discrete time designs by scheme I seem to be acceptable in the liiear 
simulations. So we expect that scheme II should be able to give more 
satisfactory compensator designs. More research effort should be speit in 
investigating them. 

In this report, we do not discuss any stochastic aspects of the 
control designs, but we believe that they are very iirportant also, ani 
a lot more should be done in this direction, especially for the discHte 
time designs which are seldom touched upon in literature. 

In developing the theoretical results for this report, we were 
just aiming at a P-I-D controller for the F-8C aircraft, and hence thi 
consideration was not made very general. We believe that one can gair 
much further and develop more adequate theory if one takes disturbanoffi 


into consideration, be they known or unknown, measurable or unmeasurarle . 
Also for simplicity and generality, of the conditions that we inaosed 


A B 
C 0 


be or a rank 


can be relaxed to simply require that the matrix 
equal to the sum of the nuitiaers of controls and outputs. Also, more 
general system or input types can be considered, e.g. , type-£ systems and 
polynomial inputs. Also a study of the robustness of the systems and its 
sensitivity to parameter variation can be very useful, as is indicated in 
Chapter V, 
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